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Abstract 

We present a nonperturbative lattice formulation of noncommutative Yang-Mills the- 
ories in arbitrary even dimension. We show that lattice regularization of a noncom- 
mutative field theory requires finite lattice volume which automatically provides both 
an ultraviolet and an infrared cutoff. We demonstrate explicitly Morita equivalence of 
commutative V(p) gauge theory with p ■ nj flavours of fundamental matter fields on 
a lattice of size L with twisted boundary conditions and noncommutative U(l) gauge 
theory with species of matter on a lattice of size p ■ L with single- valued fields. We 
discuss the relation with twisted large N reduced models and construct observables in 
noncommutative gauge theory with matter. 
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1 Introduction 



It has been suggested for some time that noncommutative geometry is a natural framework 
to describe nonperturbative string theory. This belief has been supported by the fact that 
Matrix Theory [1] or the IIB matrix model [2], which are conjectured to provide nonper- 
turbative definitions of string theories, give rise to noncommutative Yang-Mills theory on 
a toroidal compactification [3]. The particular toroidal compactification can be interpreted 
in terms of the presence of a background B field [4], which can also be understood in the 
context of open string quantization [5]. Noncommutative gauge theories possess a number of 
important properties inherent from noncommutative geometry, where there is a remarkable 
geometric equivalence relation on certain classes of noncommutative spaces known as Morita 
equivalence. In noncommutative Yang-Mills theory, this implies a duality between gauge 
theories over different noncommutative tori, for example, which relates a Yang-Mills theory 
with background magnetic flux to a gauge theory with gauge group of lower rank and no 
background flux. It allows one to interpolate continuously, through noncommutative Yang- 
Mills theories, between two ordinary Yang-Mills theories with gauge groups of different rank 
and appropriate background magnetic fluxes. Furthermore, in certain instances, there is the 
remarkable fact [6] that the non-abelian nature of a gauge group can be absorbed into the 
noncommutativity of spacetime by mapping a U(p) gauge theory with multi-valued gauge 
fields to a U(l) gauge theory with single- valued fields on a dual noncommutative torus. 

While most of the results concerning Morita equivalence are obtained at the classical 
level, there are strong indications that it persists in regularized perturbation theory. In 
particular, the one-loop divergences coincide [7] for noncommutative gauge theory and ordi- 
nary commutative Yang-Mills theory on M 4 after a proper rescaling of the coupling constant. 
Moreover, it can be shown [8, 9, 10] that noncommutative gauge theory is equivalent to all 
orders of perturbation theory to a twisted large N reduced model [11, 12]. Following this 
technique, gauge-invariant observables for noncommutative Yang-Mills theory have been 
constructed. A surprising new class of observables, which are expressed in terms of open 
Wilson loops, exists in noncommutative gauge theory [9, 13] in addition to those expressed 
via a counterpart of the standard closed Wilson loops of the reduced model. 

Lattice gauge theories [14] are a standard tool for nonperturbative investigations of non- 
supersymmetric gauge theories. In this Letter we will apply a lattice formulation of non- 
commutative gauge theory to study its nonperturbative properties and, in particular, Morita 
equivalence. The construction is an extension of a previous work [13] in which a unified frame- 
work was developed which naturally interpolates between the Matrix theory compactification 
and the twisted large N reduced model versions of noncommutative Yang-Mills theory. It 
was shown that a finite N matrix model obtained from a constrained twisted Eguchi-Kawai 
model is equivalent to a manifestly star-gauge invariant formulation of noncommutative U(l) 
gauge theory on a lattice of finite extent. In the following we will further investigate the 
properties of noncommutative Yang-Mills theories in arbitrary even dimension, and in par- 
ticular reconsider the lattice formulation of noncommutative gauge theory from a general 
point of view, without specifying a particular representation of the noncommutative alge- 
bra. Because the lattice provides a regularization of the theory, all the results obtained are 
rigorous in this sense. We will show that in noncommutative geometry, when one discretizes 
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spacetime, consistency of the algebra requires spacetime to be compactified as well, thereby 
automatically providing an infrared cutoff. This is similar to the correspondence between 
infrared and ultraviolet in noncommutative field theories discovered [15] in perturbation the- 
ory. We will also explicitly construct a map between the fields in commutative U(p) gauge 
theory on a lattice of size L with twisted boundary conditions (representing the existence 
of a 't Hooft flux) and noncommutative U(l) gauge theory on a lattice of size p ■ L with 
periodic boundary conditions, thereby demonstrating their Morita equivalence. 

The technique introduced also allows us to study properties of noncommutative gauge 
theory with matter. We will consider the coupling to matter fields in the fundamental 
representation of the gauge group and explicitly demonstrate Morita equivalence between 
commutative U(p) gauge theory with p ■ rif different flavours of matter on a lattice of size L 
with twisted boundary conditions and noncommutative U(l) gauge theory with rif flavours 
of matter on a lattice of size p ■ L with periodic boundary conditions. In the simplest case of 
L = 1 we recover the reduced model introduced in [16]. We discuss observables represented 
by the matter fields and rewrite them in terms of Wilson loops by integrating out the matter. 
This is the analog of the original construction [14] of observables of ordinary lattice gauge 
theory from matter field averages. We will construct, along these lines, observables which 
correspond to open Wilson loops and emphasize their role in the nonperturbative dynamics 
of noncommutative gauge theory. 



2 Noncommutative lattice gauge theory 

On a noncommutative space of dimension D, the local coordinates x^ are replaced by her- 
mitian operators x^ obeying the commutation relations [x^,^] = i0^ u , where 9^ = —6 Ufl 
are (dimensionful) real-valued c-numbers. To describe a lattice discretization, we restrict the 
spacetime points to x^ e eZ, where e is the lattice spacing. The lattice momentum then has 
the periodicity — > k fl + v — which implies an identity of operators such 

as e i(k„+2n6^/e)x„ _ ^k^x^ _ Acting by e - ^" on both sides, we conclude that e 2mi ^l e = 1 
and the momentum k^ should be restricted as 

Q^K G 2eZ (2.1) 

with ^Ofa, a D x D integer- valued matrix. Note that the restriction (2.1) on the momen- 
tum simply disappears in the commutative case 6^ u = 0, and is quite characteristic of the 
noncommutative geometry. 

The discretization (2.1) of the momentum implies that the lattice is periodic, x^ ~ 
x M + S^, v — 1, . . . , D, with period matrix E given by 

M^T, uX = ^9^, (2.2) 

where M^ v is a D x D integer- valued matrix. Thus, we have found that lattice regularization 
of noncommutative field theory inevitably requires the lattice to be compact as well. In the 
continuum limit e — > 0, the infrared cutoff goes away as \. For any given periodicity T,^ u 
and noncommutativity 6^ v in the continuum, one can construct a series of lattice theories 
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satisfying the restriction (2.2) and approaching the target continuum theory in the e — > 
limit. 

Lattice fields <f>i(x) on the D dimensional noncommutative spacetime are now replaced by 
finite dimensional operators fa. Their coordinate-space representation <f>i(x) can be obtained 
using a map 

fa = Y,H*)U*) (2-3) 

X 

with 

A(x) = 1 - T — [ ft (^) m J e^o 6 ^™^ e -2 7 ri(S- 1 ) M ,m^, = ? ^ 

where the sum in (2.3) runs over all lattice points modulo the lattice periodicity, and the 
sum in (2.4) goes over all integers modulo the periodicity \ T, Ufl . The operators 

= f&ip- 1 )^** (2.5) 

satisfy the commutation relations 

Z„Z V = e" 2 ™ e - Zj^ , (2.6) 
where the dimensionless noncommutativity parameter 

9^ = 27r(£-V# Ap (E- 1 ), p (2.7) 
is necessarily rational-valued on the lattice, since the restriction (2.2) implies that 

= ^S mA a , (2.8) 

should be an integer-valued matrix. The normalization of the trace of operators is fixed by 
Tr A(x) = 1. The collection of operators A(x) form an orthonormal set and the inverse 
of the map (2.3) is given by fa(x) = Tr(faA(x)). The product fafa of two operators has 
coordinate space representation given by the lattice star-product 

Tr (fafa A(x)) = \ V fa(y) fa(z) e ^(e-^M^y,K^^) <^ * ( 2 .9) 
V / det — X] 



y,z 



where we have assumed that (M -1 )^ is an integer- valued matrix 1 . 

In order to construct a lattice formulation of noncommutative Yang-Mills theory, we need 
to maintain star-gauge invariance on the lattice. As in the case of ordinary lattice gauge 
theory [14], this is achieved by putting the U(n) gauge fields on the links of the lattice [13]. 
This determines a unitary operator 

U, = J2M*)®U^x) , (2.10) 



1 This formula is similar to Rcf. [17]. See also [18] for earlier works in this regard. 
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where U^(x) is an n x n matrix field on the lattice which is star-unitary, U^x)*!/^)^ = l n . 
The lattice action is 



X fJL^V 



[2.11) 



and it is invariant under the lattice star-gauge transformation 

U^x) h-> g(x) * ^(x) * g(x + ejj) ] , (2.12) 

where the gauge function g(x) is star-unitary, g(x)*g(xy = l n . The lattice has finite extent, 
as is already discussed, and the fields U^(x) are single-valued, U^{x + T, av 6t) = U^(x). 

3 Lattice Morita equivalence 

We will now demonstrate that the noncommutative Yang-Mills theory with action (2.11) 
is Morita equivalent to a commutative Yang-Mills theory on a lattice of different size and 
with a different gauge group, thereby generalizing previous results [8, 13] for the twisted 
Eguchi-Kawai model. We start with a commutative U(p) lattice gauge theory with 't Hooft 
fluxes. The action is 

S= ^^tT^U^U^x + ef^U^x + euyUuix)^ , (3.1) 

where U^(x) are U(p) gauge fields satisfying the twisted boundary conditions 

U^x + = Sl v (x) U^x) n u (x + e/ t ) t , (3.2) 

with period matrix E. The transition functions fl^x) are SU(p) matrices obeying the con- 
sistency condition [19] 

Q^x + E au a) fl v (x) = Z^ v Vt u (x + S aM d) fl^x) , (3.3) 

where Z^ u = Q 2m Q^/p g Z p . The antisymmetric matrix Q has elements Q^ u G Z representing 
the 't Hooft fluxes. The reason why we have chosen the Q^x) as SU(p), rather than the 
usual U(p), matrices is that the abelian magnetic flux of a U(p) gauge field, which plays a 
very important role in Morita equivalences of noncommutative gauge theories [6], arises in 
the commutative case only via the corresponding 't Hooft flux [20] . A similar consideration 
with U(p) transition matrices would lead us to the same results. 

We take the gauge choice = T M with constant T M which are called twist eaters. 

Eq. (3.3) implies the Weyl-'t Hooft commutation relations 

T^T u = e 2mQ ^l p T u T^. (3.4) 

In order to find a general solution to (3.2), we need an explicit form of the twist eaters T M 
satisfying (3.4) for SU(p) and even spacetime dimension D = 2d. For generic rank p and 
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flux matrix Q, these matrices have been constructed in Ref. [21] and have the dimension 
of the irreducible representation of (3.4) equal to p/po with integer p . The most general 
gauge field configuration U^x) satisfying the constraints (3.2) is determined by two integral 
matrices P and B, which are constructed in the Appendix, and can be represented in terms 
of the operator (2.10) as 

^ = E (il (Xf j e^<^ m ^ <g> u^m) , (3.5) 

fn \u=l / 

where u^{m) is a p x p matrix and Z'^ is given by 

D 

where £' = £P and ©' = —P L B T . Because of their dependence on the twist eating 
solutions, the operators (3.6) obey the commutation relations 

Z'^l = e- 2 ^" ■ ( 3 - 7 ) 

The two matrices T! and 0' must satisfy the general constraint (2.8), which implies that 
M = -^EB T must be an integer- valued matrix. The sum over rh in (3.5) can then be 
taken over iP modulo the periodicity ~ m M + i . 
We now introduce the map 

A ' M = IdsWi ^ fn (*,)"") e" E - 9 »""" m " e- 2 "( E "')-/"»-- (3.8) 

and decompose the operator as 

^ = E A V)®W- (3-9) 

a;' 

It follows that U'^(x') is a single- valued po x p matrix field on a periodic lattice with period 
matrix £' = £P and the dimensionless noncommutativity matrix O'. Unitarity of the oper- 
ator Up further requires that U' (x 1 ) be star-unitary. Finally, we arrive at the action (2.11) 
with reduced rank n = p and coupling constant g 2 = g 2 (p/po)- This shows that an ordi- 
nary U(p) lattice gauge theory with twisted gauge fields is equivalent to noncommutative 
U(po) lattice gauge theory with periodic gauge fields. What we have arrived at is the lattice 
analog of the well-known fact that, for p = 1, noncommutative U(l) Yang-Mills theory 
with rational-valued deformation parameters Mi , is equivalent to ordinary Yang-Mills the- 
ory with gauge group U(p) and non- vanishing 't Hooft flux. From the present point of view, 
Morita equivalence is regarded as a change of basis A (re) <-> A'(x') for the mapping between 
operators and fields. 

We will now discuss the case S = elo in the present construction. A one-site U(p) lattice 
gauge theory is just the Eguchi-Kawai model [11], and the fact that the boundary condition 
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is twisted as in (3.2) means that it is actually the twisted Eguchi-Kawai model [12]. To see 
this explicitly, we reduce the action (3.1) to a single point x = by using the constraints 
(3.2) and arrive at the action 

s = -h E z ^ tr <p> v » v i kj) (3-io) 

where = U^(0) r M , /x = 1, . . . , D, are p x p unitary matrices. This is the action of the 
twisted Eguchi-Kawai model, where the phase factor is called the "twist". Thus, the 
recent proposal [8] that the twisted large N reduced model serves as a concrete definition 
of noncommutative Yang-Mills theory can be interpreted as the simplest example of Morita 
equivalence. The possibility of such an interpretation has also been suggested in Ref. [22]. 



4 Wilson loops on the lattice 

We now turn to a description of star-gauge invariant observables on the lattice [9, 13]. We 
first define the lattice parallel transport operator U(x; C) which is specified by a contour 
C = {/ii, /i2, • • • , /in}, where /ij = ±1, ±2, • • • ± D and C/_ M (x) = U^(x — efi)^ . It is defined 
by 

U(x;C) = U^(x) -kU^(x + e/ti) * • • • -kU^ ^ + e^/X, j (4.1) 

and it transforms under the star-gauge transformation (2.12) as 

U(x\ C) ^ g(x) * U{x\ C) *■ g(x + v)^ , (4.2) 

where v — e YTj=i ftj- 

Star-gauge invariant observables are constructed out of the parallel transport operator 
U(x; C) using a star-unitary function S v (x) with the property 

S v (x) * g(x) * S v (x)i = g(x + v) (4.3) 

for arbitrary functions g(x) on the periodic lattice. The solution to this equation is given by 

S v (x) = e lk ' x (4.4) 

provided that 

= O^K + H^n,, (4.5) 

with some integer-valued vector n M . Since both k and v are discrete and periodic, there are 
the same finite number of values that k and v can take, by modding out the periodicity. 
Therefore, it makes sense to ask whether (4.5) gives a one-to-one correspondence between k 
and v. The answer is affirmative if and only if there exist D x D integer-valued matrices J 
and K which satisfy 

- SJ — 2MK = !.£), (4.6) 
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where M is the integer- valued matrix introduced in (2.2). If this condition is not met, then 
there exists v for which there is no momentum k satisfying (4.5), and for the other v there 
are more than one k satisfying (4.5). For example, let us consider the case with ±£ = Lt D . 
If L is even, then (4.6) cannot be satisfied. If L is odd and M _1 is an integer- valued matrix, 
then one can satisfy (4.6) by taking J = 1^ and K = ( - L ~ 1 - ) M~ 1 . 

By using the function (4.4), with the property (4.3), star-gauge invariant observables can 
be constructed out of (4.1) as 

0(C) = tr (n) (U(x- C)) * S v (x) . (4.7) 

x 

Its star-gauge invariance stems from (4.2) and (4.3). In Eq. (4.7), the parameter k^ in (4.4) 
can be interpreted as the total momentum of the contour C. Note that in the commutative 
case, 9^ = 0, Eq. (4.5) reproduces the fact that gauge invariant quantities are given either 
by closed Wilson loops [n v = 0) or by Polyakov loops {n v ^ 0), i.e. loops winding around 
the torus. Note also that in the commutative case, the total momentum k^ is unrestricted. 
In the noncommutative case, on the other hand, one can construct gauge invariant quantities 
associated with open loops. In particular, when the condition (4.6) is met, one can construct 
a gauge invariant quantity for each contour C on the lattice, and the total momentum should 
be specified uniquely (up to periodicity) through (4.5) depending on the separation vector 
f M of the contour. We will see in the next Section that star-gauge invariant observables 
constructed in noncommutative Yang-Mills theory reduce smoothly to ordinary Wilson loops 
in the commutative limit. 



5 Coupling to fundamental matter fields 

We will now consider noncommutative gauge theories coupled to matter fields in the funda- 
mental representation of the gauge group. One advantage of the lattice formulation is that 
it allows a hopping parameter expansion (or large mass expansion) of this theory which will 
enable us to clarify various aspects of the star-gauge invariant observables constructed out 
of noncommutative gauge fields. For simplicity, we consider a complex scalar field <p(x) as 
the matter field. Fermions can be treated in an analogous way. The action for the matter 
field is 

Scatter = ~K I ^ $ (x) * U^x) * (f)(x + efi) + C.C. \ + 1 ' (x) <j>(x) , (5.1) 

L X,/J, J X 

and it is invariant under the star-gauge transformation 

<f>(x) i-> g{x)-k<j){x) ; $(x)\-+4> l (x)*g\x) ; U^x) i-> g{x) * U^x) *g\x + efi) . (5.2) 

By integrating over the field 4>(x), we perform an expansion in the hopping parameter k. 

The effective action r e ff [U] for the gauge field U^x) induced by the integration over 0(x) 
is given by 

r efr [[/] = -lnf^ /(j^^W + U^ + ^x + eft + cc^ \ . (5.3) 



n=0 \ \ x,ft 



K = 
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Throughout this section, (• • •) refers to the vacuum expectation value for fixed gauge back- 
ground, namely we integrate over the matter fields only. Using Wick's theorem, we obtain 2 

r«[U] =Y.JTc) Y, i ^r l) (u{x-C)) , (5.4) 

C x 

where J2c denotes the sum over all closed loops on the lattice and L(C) denotes the length 
of the loop C in lattice units. In this way, we encounter observables associated with closed 
loops 3 . 

Let us now consider star-gauge invariant observables involving matter fields such as 

G[f] = ^yOrWOr)* f(x)^J , (5.5) 

which is star-gauge invariant for arbitrary functions f(x) on the lattice which can be regarded 
as the wavefunction of the composite operator (ft(x)*(f)(x). We first make a Fourier transform 
of f(x) and express G[f] as 

u \ x ' 

where t> M = 0^ u k u . We have used the fact that e lk ' x acts as a translation operator in the 
noncommutative field theory as follows from Eqs. (4.3) and (4.4). Integrating over the 
matter fields using the hopping parameter expansion, we obtain 

G[f] = J2 /» E kL(c) E tr w ( w fo °) * elk ' x ) ( 5 - 7 ) 

l C x 

where J2c now denotes the sum over all lattice loops beginning at the origin and ending at 
Vp. Thus we find the observables encountered in the previous Section which are associated 
with open loops. In the commutative case, the separation vector of the loops should 
vanish independently of the momentum k^. Thus the summation over the loop C contains 
only closed loops. The summation over k can then be done explicitly reproducing the 
wavefunction f(x). In the noncommutative case, the summation over the loop C depends 
on k,j, and we are not allowed to reverse the order of the summations. The larger k^ is, the 
larger becomes the separation vector of the two ends of the loop C. This is a characteristic 
phenomenon in noncommutative field theories. If one would like to have a higher resolution 
in one direction, say in the \x = 1 direction, by increasing ki, then the object will extend in 
the other directions proportionally to 9^ik\. On the other hand, if one considers the case 

2 A calculation of the hopping parameter expansion with fundamental fermions has been performed by 
A. Zubkov. 

3 We can also introduce n/ flavours of the matter field and take the limit n/ — > oo with k ~ n^ 1 ^ 4 . In 
that case, only the simplest closed loops of lengths 4, i.e. the boundaries of lattice plaquettes, remain in 
(5.4). The single-plaquette action (2.11) can thus be induced by fundamental matter fields, similarly to the 
commutative case [23]. 
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when f(k) has support only for finite momentum then one can take the 9^ — > limit 
smoothly, reproducing the commutative case. 

The considerations in this Section show that the star-gauge invariant observables indeed 
play the same fundamental role as ordinary Wilson loops do in commutative gauge theories. 
As the knowledge of all the Wilson loop correlators would give all the information of the 
gauge theory including matter fields, so do the star-gauge invariant observables in the non- 
commutative case. We have also seen explicitly how these observables reduce smoothly to 
ordinary Wilson loops in the commutative limit. 

6 Morita equivalence with fundamental matter fields 

In this final Section, we discuss Morita equivalence of noncommutative Yang-Mills theories 
coupled to fundamental matter fields. We show, in particular, that the twisted Eguchi-Kawai 
model with fundamental matter as constructed in Ref. [16] is equivalent to noncommutative 
Yang-Mills theory with fundamental matter fields. We use the setup of Section 3. We con- 
sider commutative U(p) gauge theory and introduce Nf flavours of matter in the fundamental 
representation. We represent it as <&(x)ij, where % — 1, . . . ,p represents the color index and 
j — 1, . . . , Nf represents the flavour index. The spacetime is discretized as x^ e eZ. For the 
present purpose, it is essential to take Nf to be an integral multiple of p. In what follows, 
we assume that Nf = p for simplicity. Then, 3>(x)jj is a p x p general complex matrix. 
The action for the matter part of the U(p) gauge theory is given by 

^matter K 

U>C ) x 

which is invariant under the gauge transformation U^x) \— > g(x)U^(x)g(x + efi) and <&(#) i— > 
g(x)<&(x). The action (6.1) has also a global U(A^) flavour symmetry i— > <&(#) g', where 
g' E \J(Nf). We impose a constraint on <&(#) and U^x) given by the twisted boundary 
conditions 

U II (x + i: au a) = T u U II (x)Tl , §(x + Y, av a) = T v $(x)Tt , (6.2) 

where £ is the period matrix. Since the transition matrix lives in SU(p), this implies a 
nonvanishing 't Hooft flux of the gauge field, as we discussed in Section 3. Notice also 
that in (6.2), the T v on the left of <&(x) represents a (global) gauge transformation, whereas 
the on the right of $(x) represents a rotation in the flavour space. This is the trick 
to introducing the fundamental representation in Morita equivalence. We have made use 
of the global flavour symmetry SXJ(Nf) to mimic the boundary condition for the adjoint 
representation. Such an idea first appeared in Ref. [24] in the context of supersymmetric 
theories. 

Solving the constraint as we did in Section 3, we find that the resulting theory is noncom- 
mutative U(po) lattice gauge theory coupled to p flavours of matter fields in the fundamen- 
tal representation, with periodic boundary conditions. We can set S = eTLo to obtain the 
twisted Eguchi-Kawai model with matter fields. This is the model introduced in Ref. [16] 
as a model which reproduces large iV gauge theory with Nf flavours in the Veneziano limit 
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Nf ~ N — > oo. Now, we have found that the same model with Nf = nfN can be interpreted 
as noncommutative U(p ) gauge theory with periodic boundary conditions including n/p 
flavours of matter fields in the fundamental representation. The one-loop beta-function for 
po — 1 and fermionic matter has been calculated in Ref. [25]. The result agrees with SU(iV) 
commutative Yang-Mills theory with Nf = nfN flavours as it should due to Morita equiva- 
lence. In the present case, there is no extra infrared singularity associated with non-planar 
diagrams [15] since one- loop diagrams with fundamental matter are always planar. 
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Appendix A General representation of twist eaters 

The explicit form of the twist eaters satisfying (3.4) for SU(p) with generic rank p, flux 
matrix Q and even dimension D = 2d can be constructed as follows [21]. Using the discrete 
symmetry of the D-dimensional torus under the SL(D, Z) geometrical automorphism group, 
Q can be represented in a canonical skew-diagonal form 



Q 



9i 



V 



-Qi 






Qd 



\ 



o I 



(A.l) 



Given the d independent fluxes g« G Z, we introduce the integers pi = gcd(qi,p), pi = p/pi, 
and qi = qi/pi- By construction, pi and q~i are co-prime. A necessary and sufficient condition 
for the existence of solutions to (3.4) is that the integer pi ■ ■ -p^, the dimension of irreducible 
representation of the Weyl-'t Hooft algebra, divides the rank p [21]. In that case we write 



p 



(A.2) 



and the twist eating solutions may then be given on the subgroup SU(pi) <S> • • • <8> SU(p^) 
SU(po) of SU(p) such that Ti_i, Tj are constructed from the Weyl-'t Hooft matrices on SU(p. 
The subgroup of GL(p, C) consisting of matrices which commute with the twist eaters T M 
then GL(p ,C). 

The general solution to Eq. (3.2) for U^x) is determined by two D x D integral matrices 



is 
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which, in the basis (A.l) where Q is skew-diagonal, read 



(Pi 




V 





Pi 



B 



Pd 





Pd) 



\ 



-bi 





b d 



\ 



/ 



(A.3) 



with integral bi obeying a(Pi + b^ — 1, i — 1, . . . , d, for some integer aj. Defining matrices 
A and Q using and (ji as with P and P in (A.3), respectively, the previous condition can 
be written as AP + BQ = TL D , which is invariant under the SL(D, Z) transformation 



A 
B 



A T AA 
A T B (A- 1 ) 



P 

Q 



A- 1 P(A- 1 ) T 

a t 0(a- 1 ) t , 



(A.4) 



where A e SL (D, Z). The flux matrix Q can be written in terms of P and Q in an SL(P, Z) 
covariant way as Q = p Q P^ 1 , and it transforms as Q i— > A T Q A. One can use this symmetry 
to rotate Q back to general form and thus find the corresponding P and B which are needed 
to construct the general solution to Eq. (3.2). 
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